On the interplay between star formation and feedback in galaxy formation simulations

We investigate the star formation-feedback cycle in cosmological galaxy formation simulations, focusing on progenitors of Milky Way (MW)-sized galaxies. We find that in order to reproduce key properties of the MW progenitors, such as semi-empirically derived star formation histories and the shape of rotation curves, our implementation of star formation and stellar feedback requires 1) a combination of local early momentum feedback via radiation pressure and stellar winds and subsequent efficient supernovae feedback, and 2) efficacy of feedback that results in self-regulation of the global star formation rate on kiloparsec scales. We show that such feedback-driven self-regulation is achieved globally for a local star formation efficiency per free fall time of $\epsilon_{\rm ff}\approx 10\%$. Although this value is larger that the $\epsilon_{\rm ff}\sim 1\%$ value usually inferred from the Kennicutt-Schmidt (KS) relation, we show that it is consistent with direct observational estimates of $\epsilon_{\rm ff}$ in molecular clouds. Moreover, we show that simulations with local efficiency of $\epsilon_{\rm ff}\approx 10\%$ reproduce the global observed KS relation. Such simulations also reproduce the cosmic star formation history of the Milky Way sized galaxies and satisfy a number of other observational constraints. Conversely, we find that simulations that a priori assume an inefficient mode of star formation, instead of achieving it via stellar feedback regulation, fail to produce sufficiently vigorous outflows and do not reproduce observations. This illustrates the importance of understanding the complex interplay between star formation and feedback and the detailed processes that contribute to the feedback-regulated formation of galaxies.Comment: 20 pages, 13 figures, accepted for publication in Ap

The origin of the Milky Way globular clusters

We present a cosmological zoom-in simulation of a Milky Way-like galaxy used to explore the formation and evolution of star clusters. We investigate in particular the origin of the bimodality observed in the colour and metallicity of globular clusters, and the environmental evolution through cosmic times in the form of tidal tensors. Our results self-consistently confirm previous findings that the blue, metal-poor clusters form in satellite galaxies which are accreted onto the Milky Way, while the red, metal-rich clusters form mostly in situ or, to a lower extent in massive, self-enriched galaxies merging with the Milky Way. By monitoring the tidal fields these populations experience, we find that clusters formed in situ (generally centrally concentrated) feel significantly stronger tides than the accreted ones, both in the present-day, and when averaged over their entire life. Furthermore, we note that the tidal field experienced by Milky Way clusters is significantly weaker in the past than at present-day, confirming that it is unlikely that a power-law cluster initial mass function like that of young massive clusters, is transformed into the observed peaked distribution in the Milky Way with relaxation-driven evaporation in a tidal field.Comment: MNRAS accepte

Observing the circumgalactic medium of simulated galaxies through synthetic absorption spectra

We explore the multiphase structure of the circumgalactic medium (CGM) probed by synthetic spectra through a cosmological zoom-in galaxy formation simulation. We employ a Bayesian method for modelling a combination of absorption lines to derive physical properties of absorbers with a formal treatment of detections, including saturated systems, and non-detections in a uniform manner. We find that in the lines of sight passing through localized density structures, absorption lines of low, intermediate and high ions are present in the spectrum and overlap in velocity space. Low, intermediate and high ions can be combined to derive the mass-weighted properties of a density-varying peak, although the ions are not co-spatial within the structure. By contrast, lines of sight that go through the hot halo only exhibit detectable HI and high ions. In such lines of sight, the absorption lines are typically broad due to the complex velocity fields across the entire halo. We show that the derived gas density, temperature, and metallicity match closely the corresponding HI mass-weighted averages along the LOS. We also show that when the data quality allows, our Bayesian technique allows one to recover the underlying physical properties of LOS by incorporating both detections and non-detections. It is especially useful to include non-detections, of species such as NV or NeVIII, when the number of detections of strong absorbers, such as HI and OVI, is smaller than the number of model parameters (density, temperature, and metallicity).Comment: Accepted for publication in MNRA

Supernovae feedback propagation: the role of turbulence

Modelling the propagation of supernova (SN) bubbles, in terms of energy, momentum and spatial extent, is critical for simulations of galaxy evolution which do not capture these scales. To date, small scale models of SN feedback predict that the evolution of above-mentioned quantities can be solely parameterised by average quantities of the surrounding gas, such as density. However, most of these studies neglect the turbulent motions of this medium. In this paper, we study the propagation and evolution of SNe in turbulent environments. We confirm that the time evolution of injected energy and momentum can be characterised by the average density. However, the details of the density structure of the interstellar medium play a crucial role in the spatial extent of the bubble, even at a given average density. We demonstrate that spherically symmetric models of SN bubbles do not model well their spatial extent, and therefore cannot not be used to design sub-grid models of SNe feedback at galactic and cosmological scales.Comment: Accepted by MNRA

Characterizing gravitational instability in turbulent multi-component galactic discs

Gravitational instabilities play an important role in galaxy evolution and in shaping the interstellar medium (ISM). The ISM is observed to be highly turbulent, meaning that observables like the gas surface density and velocity dispersion depend on the size of the region over which they are measured. In this work we investigate, using simulations of Milky Way-like disc galaxies with a resolution of $\sim 9$ pc, the nature of turbulence in the ISM and how this affects the gravitational stability of galaxies. By accounting for the measured average turbulent scalings of the density and velocity fields in the stability analysis, we can more robustly characterize the average level of stability of the galaxies as a function of scale, and in a straightforward manner identify scales prone to fragmentation. Furthermore, we find that the stability of a disc with feedback-driven turbulence can be well described by a "Toomre-like" $Q$ stability criterion on all scales, whereas the classical $Q$ can formally lose its meaning on small scales if violent disc instabilities occur in models lacking pressure support from stellar feedback.Comment: 11 pages, 5 figures, submitted to MNRA

Galaxies that Shine: radiation-hydrodynamical simulations of disk galaxies

Radiation feedback is typically implemented using subgrid recipes in hydrodynamical simulations of galaxies. Very little work has so far been performed using radiation-hydrodynamics (RHD), and there is no consensus on the importance of radiation feedback in galaxy evolution. We present RHD simulations of isolated galaxy disks of different masses with a resolution of 18 pc. Besides accounting for supernova feedback, our simulations are the first galaxy-scale simulations to include RHD treatments of photo-ionisation heating and radiation pressure, from both direct optical/UV radiation and multi-scattered, re-processed infrared (IR) radiation. Photo-heating smooths and thickens the disks and suppresses star formation about as much as the inclusion of ("thermal dump") supernova feedback does. These effects decrease with galaxy mass and are mainly due to the prevention of the formation of dense clouds, as opposed to their destruction. Radiation pressure, whether from direct or IR radiation, has little effect, but for the IR radiation we show that its impact is limited by our inability to resolve the high optical depths for which multi-scattering becomes important. While artificially boosting the IR optical depths does reduce the star formation, it does so by smoothing the gas rather than by generating stronger outflows. We conclude that although higher-resolution simulations, and potentially also different supernova implementations, are needed for confirmation, our findings suggest that radiation feedback is more gentle and less effective than is often assumed in subgrid prescriptions.Comment: 28 pages, 26 figures, accepted for publication in MNRAS. Revised to match published versio

Reply to Melott's Comment on ``Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View'' by Romeo et al

Melott has made pioneering studies of the effects of particle discreteness in N-body simulations, a fundamental point that needs careful thought and analysis since all such simulations suffer from numerical noise arising from the use of finite-mass particles. Melott (arXiv:0804.0589) claims that the conclusions of our paper (arXiv:0804.0294) are essentially equivalent to those of his earlier work. Melott is wrong: he has jumped onto one of our conclusions and interpreted that in his own way. Here we point out the whys and the wherefores

Multiple populations in globular clusters: the distinct kinematic imprints of different formation scenarios

Several scenarios have been proposed to explain the presence of multiple stellar populations in globular clusters. Many of them invoke multiple generations of stars to explain the observed chemical abundance anomalies, but it has also been suggested that self-enrichment could occur via accretion of ejecta from massive stars onto the circumstellar disc of low-mass pre-main sequence stars. These scenarios imply different initial conditions for the kinematics of the various stellar populations. Given some net angular momentum initially, models for which a second generation forms from gas that collects in a cooling flow into the core of the cluster predict an initially larger rotational amplitude for the polluted stars compared to the pristine stars. This is opposite to what is expected from the accretion model, where the polluted stars are the ones crossing the core and are on preferentially radial (low-angular momentum) orbits, such that their rotational amplitude is lower. Here we present the results of a suite of $N$-body simulations with initial conditions chosen to capture the distinct kinematic properties of these pollution scenarios. We show that initial differences in the kinematics of polluted and pristine stars can survive to the present epoch in the outer parts of a large fraction of Galactic globular clusters. The differential rotation of pristine and polluted stars is identified as a unique kinematic signature that could allow us to distinguish between various scenarios, while other kinematic imprints are generally very similar from one scenario to the other.Comment: 22 pages, 16 figures + appendix. Accepted for publication in MNRA

A Systematic Look at the Effects of Radiative Feedback on Disc Galaxy Formation

Galaxy formation models and simulations rely on various feedback mechanisms to reproduce the observed baryonic scaling relations and galaxy morphologies. Although dwarf galaxy and giant elliptical properties can be explained using feedback from supernova and active galactic nuclei, Milky Way-sized galaxies still represent a challenge to current theories of galaxy formation. In this paper, we explore the possible role of feedback from stellar radiation in regulating the main properties of disk galaxies such as our own Milky Way. We have performed a suite of cosmological simulations of the same $\sim10^{12} {\rm M}_{\odot}$ halo selected based on its rather typical mass accretion history. We have implemented radiative feedback from young stars using a crude model of radiative transfer for ultraviolet (UV) and infrared (IR) radiation. However, the model is realistic enough such that the dust opacity plays a direct role in regulating the efficiency of our feedback mechanism. We have explored various models for the dust opacity, assuming different constant dust temperatures, as well as a varying dust temperature model. We find that while strong radiative feedback appears as a viable mechanism to regulate the stellar mass fraction in massive galaxies, it also prevents the formation of discs with reasonable morphologies. In models with strong stellar radiation feedback, stellar discs are systematically too thick while the gas disc morphology is completely destroyed due to the efficient mixing between the feedback-affected gas and its surroundings. At the resolution of our simulation suite, we find it impossible to preserve spiral disc morphology while at the same time expelling enough baryons to satisfy the abundance matching constraints.Comment: accepted to MNRA